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An absolute measurement of the viscosity of di(2-ethylhexyl) sebaeate with a torsional pendu- 
lum viscometer. The apparatus and calibration procedures are discussed and a detailed study of 
errors is made. The systematic error was estimated to be less than ±0.07 percent and the variability 
was shown to be less than ±0.07 percent. There is an unexplained difference of 0.3 to 0.4 percent 
between these measurements and measurements with a capillary instrument. 
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1. Introduction 

In calculating the accuracy of any physical measure- 
ment, it is not possible to be sure that all systematic 
errors have been taken into account. Moreover, the 
more the precision is improved, the greater is the pos- 
sibility that an overlooked error is significant. A 
comparison of measurements made by different 
methods, each of which insofar as possible gives rise 
to quite different errors, is the most direct evidence 
of the completeness of the error estimates. For this 
reason, some years ago we attempted a measurement 
of viscosity by an absolute method and with errors 
not exceeding 0.1 percent in order to compare with 
measurements of presumed equal accuracy made with 
a capillary viscometer. 

We chose to make the measurement with a viscom- 
eter consisting of a spherical shell suspended from a 
torsion wire and filled with the fluid whose viscosity 
is to be measured. The sphere is excited in torsional 
oscillation and the subsequent motion observed. In 
principle the viscosity of a Newtonian fluid may be ob- 
tained from either the period or the damping of this mo- 
tion. In 1959 [l] 1 an analysis of the absolute torsional 
pendulum viscometer was published giving design 
criteria for the sphere and wire suspension which 
made the period sensitive to viscosity. That paper 
pointed out that modern techniques of measuring 
time intervals with great accuracy make this a promis- 
ing method of measuring viscosity. 

This paper describes the apparatus and the pro- 
cedures used in calibrations and measurements. 
Difficulties encountered in obtaining a high degree of 
accuracy are discussed. The most reliable measure- 
ments reported are those obtained for Octoil S, a 
fluid with nominal kinematic viscosity of 19 centi- 



* Retired from the Institute for Basic Standards. National Bureau of Standards. Wash- 
ington. D.C. 20234. 

1 Figures in brackets indicate the literature references at the end of this paper. 



stokes (cSt). These measurements were subsequently 
compared with an absolute measurement with a capil- 
lary instrument [2]. Measurements of a few other 
fluids were also reported with error estimates, pri- 
marily as a consistency check. A few measurements 
by the method of decrement are reported, but without 
complete error estimates. 

2. Apparatus and Method of Measurement 

The apparatus was designed to exhibit relatively 
high dependence of period on viscosity for a range 
from 10 to 150 cSt. Above and below this range, the 
dependence is markedly weaker. Figure 1 is a photo- 
graph of the complete apparatus showing the hollow 
sphere hanging above the thermostatted and vacuum 
tight housing within which it normally operates 
and the optical devices and timer used in measuring 
the period. Figure 2 is a schematic cross section of 
the thermostatted housing showing the sphere hang- 
ing on the torsion wire in the normal operating posi- 
tion. In operation, the sphere is excited in torsional 
oscillation, and the period is measured. The kine- 
matic viscosity is then given by: 
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where v is the kinematic viscosity, K is the torsion 
constant of the suspension, R is the internal radius of 
the hollow sphere, T is the period of the oscillation, / 
is the moment of inertia of the empty sphere, and /<> is 
the moment of inertia of the test fluid considered as a 
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FIGURE 1. Photograph of the torsional pendulum viscometer. 

The sphere has heen withdrawn above the thermal chain her. 

rigid sphere. The dimensionless quantity tj is a 
mathematical function of the two dimensionless 
variables r and a. An accurate determination of vis- 
cosity requires measurements of all the physical 
quantities entering into these equations. The cal- 
culation of i] (r, a) can be done to the necessary 
accuracy by numerical methods using a modern digital 
computer. Details of this calculation are given in 
[1], (where J] is used for our r)). 

2.1 . The Torsion Constant of the Suspension 

The torsion wire was of an alloy of 92 percent 
platinum and 8 percent tungsten. This alloy was 
shown by Kestin, et al. [3] to be especially suitable 
for use as a torsion wire by reason of its low internal 
friction and stable elastic constant. The wire was a 
51-cm length of 0.051-cm diam which was heat treated 
(to about 1300 °F) to relieve stresses after which the 
ends were hard soldered into threaded fittings. 

The calibration of this viscometer requires an ex- 
tremely accurate evaluation of the torsion constant 
K. It was necessary, therefore, to measure the con- 
stant at various temperatures and for various tensile 
loads on the wire. For this purpose, a number of 
cylinders of calculable moment of inertia were made. 
The torsion constant of this wire was measured by 
timing the oscillations of these cylinders. Steel and 
brass stock were found to be unsuitable 2 for making 
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2 Several measurements of the torsion constant were made with a cylinder machined 
down in radius between each measurement. When plotted against the weights of the 
cylinders, the measurements for steel and brass cylinders did not fall on the same smooth 
curve as did measurements using aluminum or copper cylinders. 



FIGURE 2. Cross-section of thermal chamber of viscometer with 
sphere in operating position. 

these rods, presumably because the material was not 
sufficiently homogeneous. Aluminum alloy 6061 and 
oxygen-free high-conductivity copper were found to 
be suitable materials. Six cylinders were used as 
shown in table 1, with moments of inertia ranging 
from about 11,000 to L200 g • cm 2 and weights ranging 
from 2,100 to 510 g. 

The measurements of period were made with the 
cylinders oscillating in torsion while suspended in 
the position of the sphere shown in figure 2. The air 
pressure was reduced to less than 0.2 mm of mer- 
cury within the thermal jacket to eliminate any sig- 
nificant effects of the atmosphere on the period (see 
section 4). The temperature of the torsion wire was 
taken to be that indicated by the mercury in glass 
thermometers through the thermal jacket shown in 
figure 2. The temperature was controlled for a day or 
more before measurements were made to insure that 
the apparatus was in thermal equilibrium. The period 
measurements were made with the optical and timing 
systems used in the viscosity measurements. These 
systems are described in section (2.4). 
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TABLE 1. (Calibrating rods 
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inches 
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grams 


g • cm 2 


1 


6061 alloy 


1.717794 


511.31 


513.43 


1221.80 


2 


6061 alloy 


1.758891 


535.58 


537.55 


1341.14 


3 


6061 alloy 


3.010187 


1574.36 


1576.52 


11520.30 


3a 


6061 alloy 


2.655604 


1224.80 


1226.96 


6978.06 


4 


OFHC Copper 


1.758936 


1773.51 


1780.19 


4441.65 


5 


OFHC Copper 


1.907257 


2093.75 


2099.75 


6159.75 



The torsion constant of the wire was found to con- 
form to the equation 

K= 11471.5-0.0208 #-2.6(^-25) 

where K is the torsion constant in dyne cm, g is the 
tensile load on the wire in grams and t is the tempera- 
ture in degrees Celsius. 

This expression for the torsion constant is accurate 
to better than ±0.015 percent through the necessary 
range of temperature around 25 °C and through the 
range of tensile load spanned by the cylinders in table 
1. The dimensions of the cylinders were measured at 
20 °C to ±20 microinches and the masses to ±10 
milligrams. Even for measurements made with the 
smallest cylinder we estimate that the errors in these 
measurements contribute an error of at most ±0.004 
percent of the torsion constant. We estimate a further 
error of ± 0.002 percent from a possible error of, at 
most, 10 percent in the calculation of the temperature 
correction of the moment of inertia. The largest error 
in the calculation of the moment of inertia comes, 
however, from the calculation of the correction for 
the drilled holes, brass fittings, wax used to affix the 
mirror, etc. This correction amounts to 1.3 g • cm 2 . 
A 10-percent error in this calculation, which is a reason- 
able but generous estimate, would cause an error of 
±0.010 percent in the moment of inertia of the smallest 
cylinder. Errors of less than 5 /xs in the measured 
period and errors of less than 0.05 °C in the tempera- 
ture assigned to the wire caused a negligible error in 
the torsion constant. 

2.2. The Construction of the Sphere and the 
Determination of Its Radius 

Initially, a sphere was constructed of aluminum by 
machining two hemispherical cavities threaded to 
screw together at the equator. There was no way to 
be sure, however, that the fit at the threads was exact 
enough to prevent the flow of a significant amount of 
fluid into the seam between the hemispheres. Ulti- 
mately a nickel sphere was made without this seam 
by electroplating over a spun aluminum form. The 
form was machined to a sphere and ground to a tol- 
erence of 0.001 in, then fitted with nickel polar hubs 
concentric to the same tolerance. It was then plated 
with a flash plating of copper followed by a plating 
of nickel to a thickness of about 0.1 in. The aluminum 



form was removed by washing with caustic solution 
flowed through holes in the polar hubs. The exterior of 
the sphere was then turned down so that the final 
sphere had a wall thickness of 0.085 in and a nominal 
diameter of 10 cm. A rough calculation shows that 
these walls are easily stiff enough to prevent a signif- 
icant change of shape under load. 

The plating process was carried out at an elevated 
temperature (45 °C to 60 °C) so that the actual internal 
radius of the sphere was larger than the radius of the 
form at 25 °C which was 4.9535 ± .0025 cm. Since the 
difference between thermal expansion coefficients 
of aluminum and nickel is roughly 1X10 -5 / °C, the 
radius of the sphere at 25 °C should be about 4.955 
cm. The calculation of the radius in this way is not 
very accurate, however, and a more direct measure- 
ment is needed. 

For this purpose, the mass of distilled and boiled 
(approximately air-free) water needed to fill the sphere 
at 25 °C was measured by weighing the sphere empty 
and full. The level of the water was adjusted to a 
fiducial mark in the sight glass at the top of the sphere 
while it was immersed in a bath of oil thermostatted at 
25.00 °C. A value of 0.9970478 g/cm : < was used for the 
density of water and account was taken of air buoyancy 
and the volume of water in the sight glass and in the 
holes in the polar hubs. The measurement was re- 
peated with two oils of low vapor pressure using 
values of density as measured by the Metrology Sec- 
tion of NBS. These measurements give a volume of 
the sphere at 25 °C of 509.515 cm { ± 0.002 percent, 
which is the volume of a sphere of radius 4.9548 cm. 
The uncertainty of ±0.002 percent assigned includes 
the spread of a number of weighings. The uncertainty 
in the density of water is negligible. 

A value for the equatorial radius of the sphere can 
be found by filling the sphere with an oil of high vis- 
cosity and measuring the period. This is approximately 
equivalent to measuring the moment of inertia of a 
rigid sphere. Oil N with a viscosity of about 10 St was 
used for this purpose. The moment of inertia was 
calculated from eq (1) using the viscosity meas- 
ured by a capillary instrument. This value of viscosity 
need not be extremely accurate since a one percent 
error will lead to an error in effective radius of only 
about 0.005 percent. ;{ The equatorial radius measured 
in this way was 4.95525 cm. The difference between 
this and the radius of a sphere equal in volume to that 
measured is less than 0.01 percent and this difference 
may be taken as the uncertainty in the radius of the 
sphere. In fact, the evidence is that the figure of the 
"sphere" was close to that of a prolate spheroid with 
the equatorial radius slightly greater than the polar 
radius. Consequently, in all subsequent calculations 
we used a moment of inertia for the fluid calculated 
from the mass of the measured volume of fluid and the 
measured equatorial radius. 

2.3. The Moment of Inertia of the Empty Sphere 

The moment of inertia of the empty sphere was found 



:i In fact, the independent value of the viscosity of N oil measured by the decrement 
method (see section 4) was within a half percent of the capillary value. 
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to be 4074.45 g- cm 2 ±0.022 percent by timing the 
period of torsional oscillations at 25 °C in vacuum. 
This value is based on the average of a number of 
observations of the period of the empty sphere in a 
vacuum and on the previously determined value of K. 
An error of ±0.007 percent corresponds to the spread 
of a number of measurements made at different times. 
It is accounted for, in part, by the condensed or 
absorbed material on the sphere wall which apparently 
varied each time the sphere was subjected to a vacuum 
(the measured weight varied by about ±0.005 percent), 
and, further, by small changes in the mounting of the 
mirror. There is also a possible error of 0.015 percent 
associated with the value of the spring constant K. 
This error will not contribute to the scatter in the obser- 
vations of period. The actual timing errors involved in 
measuring the period of the empty sphere were com- 
pletely negligible, since they amounted to only 2 X 10 4 
percent of the period (see section 2.4). 

Equation (1) applies to a sphere of fluid contained in 
a rigid shell of moment of inertia /. Thus, in calculating 
/, the effective moment of inertia of the small amounts 
of fluid in hub and sight glass should be added to the 
moment of inertia of the empty sphere. In practice, 
these corrections for the nickel sphere are negligible 
because the small cylinders of fluid are all on the axis 
of rotation. For the original aluminum sphere, however, 
the unknown amount of fluid which leaked into the 
equatorial seam could produce a sensible effect. 

2.4. The Period of the Oscillations 

Measurements of period were made using an 
electronic timer with an optical lever and photomulti- 
plier. Light from the 0.005-in diameter straight filament 
of a shielded lamp shines through a window and lens 
(Q in figure 2) to a small galvanometer mirror (P in 
figure 2) fastened to a brass fitting on the lower pole 
of the sphere. The filament image is reflected through 
the lens and window to a horizontally adjustable 
vertical slit 0.005 in wide in an opaque film on the 
outside of a slotted black box enclosing a photo- 
multiplier tube. The output from the photomultiplier 
tube is wired to the timer and to an oscilloscope. The 
photomultiplier and electronic timer arrangement was 
checked against the NBS standard frequency and 
found to measure time intervals to within 5 ^s. Errors 
from this source are completely negligible, since only 
intervals of several seconds duration were measured. 
The principal difficulties of the measurement of period 
are to be found in the mechanical arrangements. 

Actually, half periods were measured by timing the 
intervals between zero crossings of the light beam. 
At the beginning of the oscillations, there is a transient 
effect, referred to as a "zero-point drift" in the analysis 
of the apparatus [1], which dies out in a few oscilla- 
tions. This effect is illustrated in figure 3 which shows 
some data for the aluminum alloy sphere filled with 
water. (The effect is especially large with low viscosity 
fluids and spheres of low moment of inertia.) The 
time interval between zero crossings is plotted against 
the number of the zero crossing at the beginning of 
the interval. The transient effect causes the even 



half periods to be larger than the odd ones. This 
effect can be used to adjust the slit, first by rough 
mechanical adjustment and finally by setting the trigger 
voltage of the timer. When the even and odd half 
periods approach the same limiting value, as they 
do in figure 3, the slit is set to give a correct zero 
position. 

Figure 4 shows a similar plot for the nickel sphere 
filled with Octoil, but for two different starting ampli- 
tudes. No effects of starting amplitude are evident 
although the starting effect of the previous paragraph 
is present. This fact is important in establishing that 
secondary flows due to inertial effects in the fluid 
are not causing a measurable effect on the period. 
Figure 5 shows a similar plot for N oil. For this liquid, 
with a viscosity some 50 times that of Octoil, a pro- 
nounced amplitude effect can be observed. The data 
after starting light amplitudes 4 of 12 to 36 degrees 
were obtained at the same temperature and pressure 
within a forty minute period of time, and these data 
show an increase of sphere period with pass number 
until the amplitude has decreased to about 8 degrees. 
The data for a starting amplitude of the light of 8 
degrees (obtained at a different ambient pressure 
and slightly different temperature) shows no change 
of period with pass number. 5 In measuring the period 
of the sphere filled with Octoil, half periods before 
the 5th pass are not used because of excessive "start- 
ing effect" and half-periods beyond the 9th pass 
were not used because of the very small amplitude. 
Thus, the average of half periods starting with the 5th, 
6th, 7th, and 8th passes (having light beam amplitudes 
in the range from 6 down to 3 degrees at the start 
of timing), were used in the calculation of viscosity. 

Some measurements of period were made with the 
apparatus at atmospheric pressure. In such cases a 
correction to the observed period was ascribed to the 
air drag on the outer surface of the sphere. This correc- 
tion was determined from data on the period of the 
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FIGURE 3. Period from time between zero-crossings plotted against 
number of zero-crossing at beginning of pass. 

This example of water in the aluminum sphere illustrates the transient effect. 



4 The angular amplitude of the oscillations of the light beam are twice those of the sphere. 

'Note that the experimental scatter for these measurements on a liquid with viscosity 
well outside the design range of the sphere has increased from less than 0.1 percent to 
about 0.5 percent. 
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FIGURE 4. Period plotted against pass numberfor Octoil S showing 
transient effect but no amplitude effect. 

sphere while empty and while filled with air-free oil 
N at pressures ranging from 0.2 mm of Hg to atmos- 
pheric'." The data obtained with the empty sphere 
were interpreted as having 48 percent of the air drag 
correction caused by the air within the sphere, since 
that is the percentage of surface area exposed in the 
internal cavity. With these data a correction could 
be calculated for a given period and air pressure by 
linear interpolation. At 760 mm of Hg these correc- 
tions amounted to 550 ^s. for Octoil and 470 jjls. for 
Noil. 

Vibration from nearby traffic could cause pendulum 
motions of the sphere and interfere with the measure- 
ment of period. A needle from the lower polar hub 
of the sphere (or cylinder when calibrating the torsion 
wire) extending into a container of high viscosity 
mineral oil (about 10 St at 25 °C) effectively damped 
the swinging motion with no measurable effect on 
the period of the oscillation. As a further precaution, 
the whole apparatus was mounted on a base con- 
sisting of two 30-in concrete cubes (visible in figure 1) 
connected by four heavyweight steel railroad rails 
embedded in the concrete at each end. Vibration- 
insulating pads of laminated rubber and cork composi- 
tion under the corners of each concrete cube separated 
the base from the concrete floor. 

2.5. Miscellaneous 

Octoil S is a commercial di(2-ethylhexyl) sebacate 
that has been highly purified by molecular distillation. 
It was chosen because its low vapor pressure (about 
6X10~ 7 mm Hg at 25 °C) made it suitable for use 
under vacuum conditions down to about 10 : * mm Hg. 
and its viscosity (about 19 cSt at 25 °C) is within the 
10 to 150 cSt design range of the spheres. 

Oil N is a Pennsylvania bright stock similar in vis- 
cosity to the discontinued SAE 70 motor oil. It was 
chosen because its low volatility under vacuum and 
its viscosity (about 1000 cSt) made it suitable for 
determining the effective inner radius of the sphere. 
(See section 2.2) The viscosity is low enough to permit 
filling and emptying of the sphere with reasonable 
ease and is high enough so that the period of oscilla- 



tion of the filled sphere depends only slightly on 
viscosity. 

Pressures below 1 mm of mercury were measured 
by a McLeod gage. Pressures between 1 mm and 
1 atm were determined by suitable manometers, one 
of which is shown in figure 1. Atmospheric pressure 
was measured by a conventional mercury barometer. 
A mechanical vacuum pump was used to obtain the 
pressures less than atmospheric. 

Mercury-in-glass thermometers (K in figure 2), 
graduated in tenths of degrees Celsius and calibrated 
in hundredths of degrees, were used for measuring 
the temperature inside the housing at appropriate 
positions. The calibration took into account the hori- 
zontal position of the thermometers and corrected for 
the ambient pressure at the bulb. 

Whenever the sphere was filled, it was immersed 
in a constant temperature oil bath at 25.00 ± .01 °C 
while the oil level in the sphere sight glass was ad- 
justed to a fiducial line. A syringe with a long needle 
was used for adding or removing oil. The bath temper- 
ature was controlled by a large-bulb mercury-in-glass 
thermoregulator with fixed enclosed platinum points. 
The sphere and sight glass could be used as a dilatom- 
eter to calculate changes in density associated with 
temperature and pressure variations. 

When the wire constant was being determined, 
when the moment of inertia of the empty sphere was 
being measured, or when the sphere was being used 
for viscosity determinations, oil was pumped from 
the constant temperature oil bath through the copper 
coil surrounding the sphere housing. The thermoregu- 
lator in the oil bath was set at 24.7 °C to allow for a 
0.3 °C temperature rise in the oil line following the 
pump, while a similar thermoregulator set at 24.9 °C 
and enclosed in a perforated heavy-wall copper tube 
located at the base of the housing, near the copper 
coil, controlled the room air conditioner. This arrange- 
ment maintained the ambient temperature of the 
housing at 25.0 °C. 
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FIGURE 5. Period plotted against pass number for N oil showing 
amplitude effect. 

The data at 8° light amplitude was taken at temperature different than for the other data. 
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3. Calculation of 



Viscosity 
Period 



from Measured 



The equation giving kinematic viscosity (= viscosity/ 
density), in terms of the measured period, T, and 
parameters of the sphere is repeated here, with b, the 
equatorial radius, replacing R, the radius of a perfect 
sphere of the same volume (see section 2.2) 



v = b 2 \ — 7/(a,r) 

y jo 



a = f,/ = |pF6 2 

i() D 



T = 7 



Air-free Octoil S at 0.24 mm Hg, 24.97 °C. p = 0.91014 

g/cm 3 
Observed period, T= 4.297637s. 
The parameters of the sphere listed above are then: 

/ = | Vpb 2 = 4554.660 g • cm 2 

T = —J^-= 1.084750 
2tt V/o 

a =///« = 0.894566 

In the table of calculated values of r is a line with 
^ = 0.89450, r=1.084737, ^ = 0.0050100, and in this 
region of the table At5/At = 0.07042, At}/ Ac* =-0.03380. 
xTius, a linear interpolation gives us for the tj cor- 
responding to the observed r and a: 



b = equatorial radius = 4.95525 cm 8 b/b = ± 1.0 

Xl0- 4 (sec.2.2) 
K = torsion constant of wire = 11455.54 dyn-cm 

8X/£ = ±1.5xlO- 4 (sec. 2.1) 
/= measured moment of inertia of empty sphere 

= 4074.445 g • cm 2 81/1=8 KfK 
J = moment of inertia of liquid in sphere moving 

as rigid body 
V= volume of sphere excluding sight glass and 

hubs = 509.515cm 3 
p = density of test fluid 
8Iolh=8plp + 8VIV+28blb 
8V/V=±2xlO-» 
Sp/p=±2XlO" 5 

The above uncertainties refer to systematic 7 errors 
only, that is those associated with constants of the 
apparatus and the density of the test fluid. The syste- 
matic error in period was less than five microseconds 
and the periods were two seconds or longer. That in 
weighings was less than 1 mg in values of 270 grams 
or greater. Neither contribute significantly to the 
total systematic error, though there were differences 
between replicate determinations due to variations in 
the mass of the sphere (see section 2.3). 

The functional relationships between t), a, and r 
comes from the solution of the equations of motion 
assuming that the fluid motion has no component in 
planes containing the axis of rotation of the sphere 
(no secondary flow). It is an implicit relation which 
can be used to calculate the r's corresponding to given 
values of r) and a. A set of such solutions was computed 
and tabulated with a fine enough grid so that linear 
interpolation was permissible. Thus, using an observed 
period T (from which r could be calculated) and a 
known a, a corresponding value of 77 could be found 
from the computed tables by interpolation and a 
corresponding value of v obtained. 

Example of calculation of viscosity from period: 



7 These errors will not produce a scatter amonj; replicate determinations of the viscosity 
of a given fluid. 



fj = 0.0050086 

from which we calculate 

v = yb 2 (Klh) ll2 = 19.504 est. 

The systematic error associated with this value 
may be calculated in terms of the independent con- 
stants of the experiment b,K,p, and Fas 

8v_ / _^d2?_ 9 a dr}\ (8b 8p 
v \ i) dr r) da J \ b p 

+ 2\/^ dr rjda) \ K V ) 

= ±6.2X10- 4 , 

this 0.06 percent is the sum of the absolute values of 
the four independent contributions, that associated 
with K giving 0.03 percent, that with b giving 0.02 
percent, while the other two contributions are each 
less than 0.01 percent. 

In order to compare this value of viscosity to that 
obtained at atmospheric pressure with a capillary 
instrument, one must know the effect of pressure 
on the viscosity. At 25.00 °C and pressures below 10 
mm of Hg, the viscosities of Octoil S and N oil are 
lower than that at atmospheric pressure by 0.20 
percent and 0.30 percent respectively. This value can 
be calculated for Octoil from published measurements 
of viscosity and density as a function of pressure for 
both Octoil and Plexol, a less refined commercial 
grade of the same chemical. These liquids are labeled 
samples 1— A and 19-D respectively in an ASTM 
report on pressure effects on lubricating fluids [4]. 
The effect of pressure on the viscosity of N oil was 
measured with relative viscometers in our laboratory, 
and the effect on the density of both oils measured 
using the sphere as a dilatometer. 

Corrections of p and v for the effect of being air-free 
were also needed for both oil N and Octoil S. For air- 
free oil N at 25.00 °C the density was greater by 
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0.00007 g/cm 3 and the viscosity was greater by 0.038 
St (0.36 percent) than that of air saturated oil N; for 
Octoil S the corresponding density difference is 
0.000010 g/cm : * and the viscosity difference is 0.00068 
St (0.35 percent). With oil N, these corrections were 
estimated from relative measurements with a capil- 
lary viscometer and from weighings of a known volume 
in the NBS density laboratory. Air-free Octoil S (much 
lower in viscosity than oil N) apparently dissolved 
air rapidly while being transferred from the sphere 
to capillary viscometers or volumetric flasks. Accord- 
ingly, for Octoil S these corrections were estimated 
from relative measurements made with the sample 
in the sphere and at atmospheric pressure. 

4. Determination of the Viscosity of N oil by 
Measurement of Decrement 

Since the viscosity of oil N is used in the determina- 
tion of the radius of the sphere (see section 2.2) the 
apparatus was adapted to give an independent value of 
this viscosity through measurement of the decrement 
of the torsional oscillation. 

The apparatus was designed for the determination 
of viscosities by measurement of periods, and origi- 
nally, no provision was made for determining decre- 
ment. We found, however, that quite accurate decre- 
ments could be measured through timing successive 
passes of the light beam through a small interval near 
the zero rest position. The single 0.005-in slit in front 
of the photomultiplier tube was replaced with two 
similar slits with their centers 0.065 inches apart. 
At rest, the light beam triggering edge was adjusted 
midway between the two slits. In motion, the light 
beam triggering edge started the timer as it passed 
the center of one slit and stopped it as it passed the 
corresponding position in the second slit. The 0.065-in 
interval was chosen to be large enough for good ac- 
curacy of timing and small enough so that the average 
velocity over the interval would be nearly the same as 
the instantaneous velocity at the zero rest position. 

For computing the viscosity from decrement, the 
ratios of timings for the 0.065-in light beam travel 
near the zero rest position, for successive passes in 
the same direction, were used. The ratio was measured 
for beam travel in each direction and then averaged 
to minimize the error caused by any inexactitude in 
centering of the interval. If t n denotes the time for nth 
traverse of the interval, D denotes a decrement defined 
by 

D= (logetnlt n +l)lTT. 

This decrement must be corrected for the damping 
caused by the internal friction of the torsion wire and 
the viscous "air-drag." 

Decrement measurements were made with the 
empty sphere oscillating in air from atmospheric 
pressure down to 0.008 mm Hg. The damping decreased 
with decreasing pressure at pressures above 1.0 mm 
Hg. remained almost constant while the pressure was 
between 1.0 and 0.1 mm, and again decreased as the 



pressure was decreased below 0.1 mm. Measurements 
of the period of the empty sphere showed no change 
throughout the range of pressures from 1.0 to 0.01 
mm Hg. Presumably, the effects of pressure changes 
down to pressures of 1.0 mm of Hg are primarily 
through changes of the density of air while those for 
pressures below 0.1 mm of Hg are associated with the 
changing viscosity or mean free path of the air. Data 
were obtained with oil N35 and Octoil S at pressures 
down to about 0.1 mm Hg. 

In the range of 0.1 to 1 mm Hg ambient pressure, 
the magnitude of the corrections to the decrement, 
D, were estimated as: 

correction due to internal friction of wire 0.000067 
correction due to air drag on sphere 0.000031. 
These should be compared to a total decrement for 
oil N of about 0.04. 

The viscosity of oil N obtained by this decrement 
method was within 0.5 percent of the value obtained 
by capillary methods relative to the currently accepted 
value for the viscosity of water at 20 °C. 

5. Results 

The collection of measured values of kinematic 
viscosity of Octoil S is shown in table 2. These meas- 
surements were made at various times over a period 
of about two years. Values measured with the original 
machined aluminum sphere are included and are 
marked with one or two asterisks to indicate two 
conditions of smoothness of the interior surface. 8 
Since the temperature and pressure varied slightly 
for the different measurements, we compare them in 
terms of the difference between the viscosity measured 
in the sphere and that measured in a relative capillary 
viscometer, calibrated in terms of the generally accept- 
ed value for the viscosity of water, at the same temper- 
ature and pressure [2]. Measurements were also 
made on several test fluids outside the design range 
of the apparatus to see if the deviations are consistent 
with our analysis of the instrument. These fluids, water 
and two of the standard viscosity fluids formerly 
distributed by NBS, each had a vapor pressure at room 
temperature which excluded the possibility of operat- 
ing the instrument in vacuum. Consequently, all the 
determinations of viscosity on these three fluids 
involved a period corrected to account for air drag 
(see section 2.4.). The results are shown in table 3. 
Determinations of viscosity by the decrement method 

are found in table 4. 

A systematic error of ± 0.06 percent was estimated 

for the measurement of the viscosity of Octoil S as 

noted in section 3. Calculations for water. Oil D, and 

Oil L yielded values ranging from ±0.05 to ±0.07 

percent, an insignificant variation, with an average for 

the four calculated values of ±0.065 percent. We may 

therefore assign an estimated systematic error of 

±0.07 percent to the measurement of any fluid in 

this range of viscosity (1 to 90 cSt). 



8 The aluminum sphere consisted of two hemispheres which screwed together at the 
equator. The resulting seam between the hemispheres introduced an error in the measured 
viscosity which was probably as much as().] percent. 



547 



_^ 



TABLE 2. Comparison of viscosity of Octoil S {in centistokes) measured by sphere period {absolute method) and capillary {relative method 

with value calculated based on viscosity of water at 20 °C= 1.004 cSt) 





Ambient 


^sphere. 




Fluid and condition 


pressure 


from period 


— x 100 






measurements 


^cap 




mm Hg 


cSt 


percent 


Octoil S, air-saturated 


758 


19.475 


0.404 




757 


19.477 


.414 


K ap = 19.397) 


745 


19.468 


.366 




738 


19.468 


.364 




735 


19.470 


.375 


Average 


atmos 


19.472 


0.385 


Octoil S, air-free 


766 


19.541 


0.389 




766 


19.539 


.382 


(j> ca p= 19.465) 


766 


19.537 


.369 




756 


19.544 


.404 




759 


19.548 


.427 




759 


19.547 


.419 




760 


19.550 


.435 




740 


19.541 


.390 




755 


19.540 


.383 


Average 


atmos 


19.543 


0.400 


Octoil S, air-free 


0.24 


19.484 


0.297 


(j> cap = 19.426) 


.16 


19.497 


.368 




2.0 


19.497 


.363 




0.9 


19.510 


.433 


Average 


vac 


19.497 


0.365 


Octoil S, average, all 






0.388 










Octoil S,* air-free 


756 


19.533 


0.347 




756 


19.517 


.269 


K,, p = 19.465) 


756 


19.513 


.245 




755 


19.508 


.221 




753 


19.524 


.303 




753 


19.525 


.307 


Average 


atmos 


19.520 


0.282 


Octoil S,* air-free 


0.15 


19.495 


0.356 




.15 


19.488 


.321 


(*/ cap = 19.426) 


.13 


19.483 


.293 




.12 


19.492 


.342 




.11 


19.495 


.354 


Average 


vac 


19.491 


0.333 


Octoil S,* average, all 






0.305 










Octoil S,** air-free 


756 


19.570 


0.538 




756 


19.572 


.549 


(* cap = 19.465) 


760 


19.556 


.465 


Average 


atmos 


19.566 


0.517 


Octoil S,** air-free 


1.0 


19.522 


0.493 




0.22 


19.524 


.504 


(i> cap = 19.426) 


0.22 


19.528 


.524 


Average 


vac 


19.525 


0.507 


Octoil S,** average, all 






0.512 











* In aluminum alloy sphere as machined. 

**In aluminum alloy sphere after polishing interior. 
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TABLE 3. Comparison of viscosities (in centistokes) measured by sphere period (absolute method) and capillary (relative method with 

value calculated based on viscosity of water at 20 °C = 1.004 cSt). 





Ambient 
pressure 


^sphere. 

from period 


^sphere "eap^ jqq 


Fluid and condition 








measurements 


^cap 




mm Hg 


cSt 


percent 


L31, air- saturated 


752 


86.991 


0.275 




752 


87.020 


.309 


(vca P = 86.752) 


752 


86.996 


.281 




755 


87.047 


.340 




758 


86.965 


.246 




758 


86.983 


.266 


Average 


atmos 


87.000 


0.286 


Water, air-saturated 


745 


0.90020 


0.813 




746 


.89992 


.782 


(i/ cap = 0.89294) 


752 


.89928 


.710 




754 


.89857 


.631 




759 


.89855 


.628 




759 


.89828 


.598 




753 


.90087 


.888 




756 


.90199 


1.013 




756 


.90196 


1.010 




754 


.90191 


1.005 


Average 


atmos 


0.90015 


0.808 


Water, air-free 


760 


0.89518 


0.251 




760 


.89544 


.280 


(v C ap = 0.89294 


759 


.89518 


.251 




754 


.90050 


.847 




754 


.90041 


.837 




752 


.90147 


.955 




751 


.90215 


1.031 




751 


.90168 


.978 




752 


.89693 


.447 




751 


.89681 


.433 




751 


.89654 


.403 




757 


.89699 


.454 




759 


.90098 


.900 




758 


.90139 


.946 




759 


.90128 


.934 




761 


.89989 


.778 




755 


.89632 


.378 




755 


.89723 


.480 




759 


.89663 


.413 


Average 


atmos 


0.89858 


0.631 


Du, air-saturated 


752 


3.1859 


0.512 




751 


3.1856 


.502 


{v C ap = 3.1697) 


751 


3.1855 


.498 




742 


3.1856 


.502 


Average 


atmos 


3.1857 


0.504 


Du, air-saturated 


756 


3.1913 


0.511 




755 


3.1885 


.422 


(p cap = 3.1751) 


754 


3.1934 


.577 




752 


3.1932 


.569 




756 


3.1937 


.585 




748 


3.1936 


.584 


Average 


atmos 


3.1923 


0.541 


D| », average, all 






0.526 
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Table 4. Comparison of viscosities {in centistokes) measured by decrement method and capillary {relative method with value calculated 

based on viscosity of water at 20 °C= 1.004 cSt). 





Ambient 


^sphere from 




Fluid and condition 


pressure 


decrement 


p — x 100 






measurements 


^cap 




mm Hg 


cSt 


percent 


N35, air-free 


762 


1054.3 


0.603 




761 


1052.8 


.457 


(v cap = 1048.0) 


753 


1053.8 


.555 


Average 


atmos 


1053.6 


0.538 


N35, air-free 


632 


1052.3 


0.460 




524 


1050.7 


.354 




303 


1050.0 


.374 




164 


1051.2 


.540 




145 


1051.8 


.606 




59 


1050.4 


.510 




21 


1050.2 


.497 


Average 






0.477 










N35, air-free 


2.5 


1050.4 


0.528 




2.0 


1050.1 


.500 


{v C ap= 1044.9) 


0.45 


1050.5 


.540 




.25 


1051.4 


.623 


Average 


vac 


1050.6 


0.548 


N35, average, all 






0.510 








Octoil S, air-free 


762 


19.391 


-0.382 


1 


762 


19.395 


- .360 


(v cap = 19.465) 


760 


19.391 


- .382 




760 


19.452 


- .067 




760 


19.432 


- .167 




758 


19.479 


.070 




754 


19.483 


.093 


Average 


atmos 


19.432 


-0.171 


Octoil S, air-free 


645 


19.460 


0.004 




531 


19.450 


- .016 




431 


19.469 


.107 




363 


19.516 


.365 




317 


19.463 


.107 




154 


19.453 


.100 




150 


19.486 


.265 




37 


19.467 


.200 


Average 






0.142 








Octoil S, air-free 


9 


19.468 


0.216 




9 


19.446 


.102 


(v C ap= 19.426) 


8 


19.427 


.005 




3 


19.422 


- .019 




2 


19.493 


.344 




0.33 


19.441 


.075 




.29 


19.423 


- .014 


Average 


vac 


19.446 


0.101 


Octoil S, average, all 






0.029 









In addition there was observed, of course, a variabil- 
ity between replicate measurements of the same 
fluid. There are primarily two factors which are 
responsible for this. First we expect small variations 
in the actual temperature of the fluid and the torsion 
wire which we estimate to be ±0.005 °C and ±0.05 



°C respectively. This will affect the viscosity of the 
fluid being measured and the torsion constant of the 
wire. The effect on dimensions of the sphere is not 
significant. The second factor, and the major one. is 
the observed variation in the weight of the sphere 
which seems to account for the variability of 0.007 
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percent in the moment of inertia of the empty sphere 
when measured in vacuum. For measurements at 
atmospheric pressure this variability in the measured 
moment of inertia increased to 0.01 percent. Pre- 
sumably these variations in / will also occur between 
measurements of viscosity and we have calculated 
the variability in the measured viscosity by making 
this assumption. 

Measurements of viscosity by the period method 
are plotted in figure 6. The results are presented in 
terms of the percent deviation of individual measure- 
ments from the results of relative viscometry with 
capillary instruments taking the viscosity of water 
as 1.002 centipoise (1.004 centistokes) at 20 °C. 

The dashed lines are centered on the mean of the 
observed points for each liquid to indicate the estimate 
of spread of errors calculated above. The agreement 
between the observed variability and that estimated 
is as close as we could expect. For water, Octoil S, 
and Oil L the agreement is almost exact. For Oil D the 
observed spread is smaller than that estimated, 
probably because of the limited number of observa- 
tions. The calculated systematic error and variability 
for period measurements on Oil N were 3.9 percent 
and 1.1 percent respectively. The spread of a few 
observations of the period with Oil N in the sphere 
were consistent with the latter value; the dashed lines 
are continued toward this range of viscosity merely 
to indicate that we expected, and observed, greater 
variability on both sides of the design range. 

In an attempt to see if the smoothness of the inner 
surface of the sphere might have any influence on the 
observed period, the original tool marks in the alumi- 
num sphere were polished off, changing the roughness 9 
from the original 20 microinches to 5 microinches. 
Of course, a new set of constants for the sphere had to 
be determined after this polishing. Table 2 includes 
values of viscosity obtained with the aluminum sphere 
both before (*) and after (**) this polishing procedure. 
The difference of about 0.2 percent is consistent with a 
reasonable set of error estimates for the aluminum 
sphere, in which we know that there is an uncertainty 
in the determination of / due to the seam between the 
two hemispheres (see section 2.3). If the roughness of 
our spheres did influence our measurements, we would 
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FIGURE 6. The difference in percent between the measurement of 
viscosity with the sphere and the relative value by capillary 
measurement. 

The zero corresponds t<> capillary measurements based <>n the current standard of 1.002 cP 
for the viscosity of water at 20 °C. 

expect greater roughness to correspond to higher 
values of viscosity. The difference observed here was 
in the opposite direction. It is clear from figure 6 that 
all the measurements of viscosity made by timing 
the period of the sphere are mutually consistent among 
themselves. It is also clear from figure 6 that our meas- 
urements of the viscosity in the design range of the ap- 
paratus (Octoil S and Oil L) are 0.3 percent or 0.4 
percent higher than the relative measurements with a 
capillary instrument taking the viscosity of water as 
1.002 cP (1.004 cSt) at 20 °C. These errors are beyond 
our estimate of the maximum possible error. We can 
only conclude that there is an unidentified system- 
atic error in one or both of these methods of measuring 
viscosity. 
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! ' Tliese values represent the arithmetical avera^r of tin 
measured with a profilometer. 
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